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so doing we discuss the prescription of overlapping singularities in more than one variable. 
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Introduction 



In recent years there has been an increasing wealth of interest in semi-inclusive deep inelastic scattering, driven 
both by crucial breakthroughs in the QCD description of these processes yj. y. LI and also by an incipient availability 
ff) ' of data encompassing polarized, unpolarized, leading baryon, and diffractive deep inelastic phenomena^- 

From the perturbative QCD standpoint, semi-inclusive deep inelastic scattering (SIDIS) brings before theorists two 
C*~) ' novel and interesting features. On the one hand, fracture functions which, in addition to structure and fragmentation 
functions, are required for the correct description of hadrons produced in the forward direction and for the factorization 
■ of collinear singularities. On the other, non-homogeneous Altarelli Parisi evolution equations, which highlights the 
Q_i interplay of the three intervening parton densities in the scale dependence of these processes 0. 

JL , Although the main features related to SIDIS, and specifically to fracture functions, have been studied at the leading 
order (LO) in QCD (order-a s in the cross-section) 0, 0] , up to now no computations had been done up to next to 
leading order (NLO) accuracy, as it is standard in the inclusive case. In particular, there were neither explicit checks 
L* • of factorization at order-a 2 . nor indications of how relevant the non homogeneous evolution might be at NLO. 
. £^ | In LO, non standard evolution effects although non negligible, are restricted to a relatively small kinematic region, 
associated to the fragmentation configurations allowed at that order[(|. This suggests to neglect these effects in many 
phenomenological analyses of polarized SIDIS 7, Q , leading baryon production and diffractive DIS 6] , provided 
some cuts on data are introduced. In NLO the above mentioned kinematical restrictions are no longer present, which 
in principle may lead to important corrections. In any case, their phenomenological relevance needs to be assessed. 

From a theoretical point of view, the computation of the SIDIS NLO corrections, and specifically the explicit check 
of factorization of collinear singularities involve also some subtleties which need close attention. At variance with 
the totally inclusive case|lfjj. where after a convenient integration over final states the remaining singularities may be 
written as distributions in only one variable times a regular function, in the one particle inclusive case at order a 2 , 
it is necessary to keep additional variables unintegrated. Consequently one, must deal with entangled singularities in 
more than one variable, corresponding for example to three particles becoming collinear simultaneously. As usual for 
semi-inclusive processes, in order to check factorization on has to keep track of the kinematical origin or configuration 
which gives rise to the singularity, which represents a non trivial additional complication and requires a detailed 
analysis of the singularity structure characteristic of the process. 

In this paper we address the above mentioned issues restricting ourselves to processes where the initial state parton 
is a gluon. This allows to analyze and answer the main issues involved skipping for the moment, and for the sake 
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of clarity, the formidable singularity structure associated with virtual corrections to quark initiated processes, which 
will be addressed in a forthcoming publication. In doing this, we develop suitable prescription rules for dealing with 
the SIDIS singularity structure. 

As result of our approach we obtain the hitherto unknown NLO non homogeneous kernels for fracture functions 
and discuss their distinctive features such as their non-factorizable dependence upon two variables. We also verify 
explicitly the factorization of collinear singularities up to order a 2 , and give the expression for the renormalized 
fracture function in terms of the bare one. In order to asses the relevance of NLO corrections we compare the effects 
of the new evolution kernels with the already known LO corrections. 

The outline of the paper is the following: in the next section we introduce the relevant kinematics and conventions 
used and we extend the 0(a s ) results for the SIDIS cross sections as required for the later factorization of collinear 
singularities at 0{a 2 a ). In the second section we discuss the computation of amplitudes and phase space integration 
of the 0{ct 2 s ) processes. There, we introduce a suitable parameterization for the phase space of the three final state 
particles, and extend some of the results given in refs. in mm for the angular integration of the corresponding 
amplitudes. In the third section we analyze the SIDIS singularity structure at order a 2 s and give details about the 
prescription recipes required for dealing with it. In the fourth we address the issue of factorization and technicalities 
associated with the convolution of distributions in many variables, present the novel NLO kernels and discuss the 
evolution of fracture functions. In the last section we present our conclusions. 



I. KINEMATICS AND 0(a s ) RESULTS 

We begin considering the one-particle inclusive process in which a lepton of momentum I scatters off a nucleon of 
momentum P, 

l(l) + P(P) -^l'(l') + h(P h )+X, (1) 

and where in addition to the emerging lepton of momentum I' , a hadron h of momentum Ph is tagged in the final state. 
X stands for all the unobserved particles. For simplicity we consider only the exchange of one photon of momentum 
q = I' — I. In order to characterize the hadronic final state, in addition to the usual DIS variables 

Q 2 = -q 2 = -{l' X B — 7pFT~ ' V=i4, S H = (P + l) 2 , (2) 

2F ■ q P ■ I 

we introduce energy and angular variables 

E h l-cos6» h 

Vh = TTTi V ' w h = o ' 3 

E (l-x B ) 2 

where Eh and E are the energies of the produced hadron and of the incoming proton in the P + q = frame, 
respectively. Oh is the angle between the momenta of the hadron and the virtual photon in the same frame. 
The corresponding cross section, differential in the final state lepton and hadron variables, can be written as0 

da ^ f 1 du f 1 dvj f 1 ( x B \ „ ( v h \ da^ 



dxB dy dvhdwh i g « 'x B u Jv h v j Jo V u J \ v i J dxBdydvjdw 



+ E t ~%P (-,(1-^1 (l-x B )-^-S(l-w h ), (4) 

i 17 l-(l-* B )„ h " y 

where the sum is over all parton species. 

In the first term of the r.h.s of eq.J3J, d&ij represents the partonic cross section for the process I + i — > I' + j + X, 
whereas /j/p and D^/j are the usual partonic densities and fragmentation functions. The variable u is related to the 
fraction of momentum of the incoming parton £ by £ = Xb/u, while Vj and wj are the partonic analogs of Vh and Wh- 
This term, represented in fig^ 1 , describes a 'current fragmentation' process in which a final state parton j fragments 
into the final state hadron h, which is produced in the same direction than j. 

In the second term of the r.h.s of eq.Q, d&i stands for the inclusive partonic cross section initiated by parton i and 
is convoluted with the fracture functions M^ip. This term, shown in figH 3 , corresponds to a 'target fragmentation' 
process, where the initial state nucleon fragments into the final state hadron and a parton, i, which participates in 
the hard scattering. In the last case, the hadron is produced in the direction of the incoming nucleon. 
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FIG. 1: (a) Current and (b) target fragmentation processes. 



The above mentioned partonic cross sections can be calculated order by order in perturbation theory and are related 

(„) („) 

to the parton-photon squared matrix elements H „„ (i, j) and H ' (i) for the i + 7 — > j + X and i + 7 — > X processes, 
respectively: 



dxB dy xbSh \ Q 

where n runs over the number of particles in the final state. Matrix elements are averaged over initial state polariza- 
tions, summed over final state polarizations and integrated over the phase space of the unobserved particles. j( n > is 
the Jacobian coming from the phase space integration and depends upon the number of final state particles n. a em 
stands for the fine structure constant and e is the electron charge. Finally, Ym and Yl are the standard kinematic 
factors for the contributions of each photon polarization and are given by, 

l + (l-y) 2 v l + 4(l-y) + (l-y) 2 



FIG. 2: Born contribution to the cross sections 

The total inclusive cross sections are well known up to order a 2 s , and more recently there have been impressive 
efforts to go beyond the NLO jlil . The corresponding complete expressions for the singular and finite pieces up to 
order a 2 s can be found in ref . [l0| . 

For the one-particle inclusive cross section, the zeroth-order in a s comes from the diagram in figElgfving, in d = 4+e 
dimensions, the qq cross section: 



da {Q) 



qq.M 



dxB dy dv dw 

with 



= c g S(l-u)S(l-v)8(w), (7) 



^{2 + e)e 2 q Y M . (8) 

2 XBSH 
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The antiquark cross section d&qq is identical to da qq whereas all the remaining cross sections vanish. The index M 
refers to the metric terms in eq.(j5|, longitudinal contributions are absent at tree level. Notice that at this order the 
quark is always produced in the backward (w — 0) direction implying that forward hadrons (w = 1) would come 
solely from target fragmentation processes, that is, those taken into account by fracture functions. 








FIG. 3: Real and virtual contributions to the a a cross sections 



The first order corrections to the one-particle inclusive cross section are also known. Expressions for the singular 
and finite terms in dimensional regularization[l5j can be found in 0] for the unpolarized case and in Q for the 
polarized one. In order to accomplish the factorization of collinear singularities at O(al), one also needs the 0(a s ) 
cross sections up to order e, for this reason we accordingly extend here the results of 0- The corresponding diagrams 
are shown in figOJ As it is explained in ref.Q, the integration region for the cross section, eq.Q), need to be splitted 
into two regions, Bl and B2 respectively, in order to account for kinematical constraints in the phase space: 



Bl = {u G [xb, x u ], »£ [a, 1], w € [0, w r }} 
B2 = {u G [x u , 1], v G [v h , 1], W G [0, w r ]} , 



(9) 



with x u ~ xb/(xb + {1 — xs)vh). The metric terms of the unpolarized cross sections can be expressed in the following 
form: 



da 



(i) 

qq(qq),M 



(1) 

qg(qg),M 



da , — \ 

gq(gq),M 



da 



da 



(1) 

qq(qq),M 
(1) 

qg(qg),M 



da (1) 

a(J gq(gq),M 



Bl 



Bl 



Bl 



B2 



B2 



B2 



Cq C e 



C q C e 



Cq Ce 



Cq C e 



'\qq,M 



eD 



(i) 

lqq.M 



pW q (u)6(l-v)6(w)+Ci 
P£l q ( u ) S(a - v) 5(1 -w) + C[H M + efl£ |M } 



-P, 

e 

2 
e 
2 
e 



(0) 



,(«) 5(v -a)6(l-w) + ~ pM g ( u ) 6(1 ~ v) 5(w) + c[^ M 



\ (Pfl q (u) 6(1 -v) + P^ q (v) 6(1 - u)) 6(w) + C^ qM + 



eD 



(i) 

lqq.M 



c q C e \- £ P(°2 q (v) 6(1 - u) 5(w) + C (1) 

2 , , , „, , ^,(1) 



2qg,M 

C e {^P^ g (u)5(l-v)6(w)+ci gqAI 



eD 



(i) 

2qg,M 



eD 



(i) 

2gq,M 



where C e is defined by 



and 



r - as f ( Q2 



(1 - u) x B 



e/2 



/r 



r(l + e/2) 



(1 



XB) 



r(i + e) 

(1 - ^) (1 - u)x B 
v(u — x B ) 



(10) 

(11) 
(12) 

(13) 
(14) 
(15) 

(16) 
(17) 



The P}°lj are the usual LO Altarelli-Parisi kernels 0], whereas the Pj^l-i are the unsubstracted ones. Expressions 
for them and for the coefficient functions can be found in Appendix B of 0. Notice that in order to match our 
notation the results of ref.0 should be multiplied by a factor S(w r — w). The coefficients D^ are explicitly given 
in Appendix I below. Similar expressions for the longitudinal cross sections, which are finite at this order, can be 
obtained from the results in ref. Q- 

Notice that at this order, the fragmenting parton can be produced in any direction, including the forward one, but all 
singular contributions come either from the backward or from the forward direction. The former terms are factorized 
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into partonic densities and fragmentation functions, whereas the forward singularities can only be factorized in the 
redefinition of fracture functions. This factorization gives rise, as we have already mentioned, to non-homogeneous 
terms in the evolution equations of fracture functions dj. Notice also that singular terms in the forward direction are 
always accompanied by a S(v — a) factor. This is a characteristic feature of the LO results that in general, will not 
be present at the a 2 order. 



II. 0{a 2 s ) AMPLITUDES AND PHASE SPACE INTEGRATION. 

In this section we outline the computation of the order a 2 amplitudes for the one-particle inclusive cross sections. At 
this order the relevant amplitudes have either two or three final state partons, related to virtual and real contributions, 
respectively. We have computed the corresponding hadronic tensors H^v in d = 4 + e dimensions, in the Feynman 
gauge, and taking all the quarks to be massless. 

Algebraic manipulations were performed with the aid of the program Mathematica an d the package Tracer 
|l8j to perform the traces over the Dirac indices. In order to obtain the one-particle inclusive cross sections, one has 
to integrate the resulting matrix elements over the internal loop momenta and over the phase space of the unobserved 
particles in the final state, which is one of the hardest and most delicate parts of the calculation. 









FIG. 4: Real contributions to the a. cross sections 



As we mentioned, in the present paper we restrict ourselves to gluon initiated processes, the corresponding 0{a 2 s ) 
real contributions (gg and gq processes) are shown in fig0] For the first of these processes, the phase space integration 
is over the momenta of the quark-antiquark pair, whereas for the second it is over the gluon-antiquark momenta. To 
perform this integrals, we choose to work in the center of mass frame of the two unobserved partons and get for the 
phase space: 



dPS 3 = -L 



f 



Q 2 T(l +e 
x(l — w) e / 2 dvdw sin 



1 - x B 



l+e/2 



X B 

1+£ fa sin e (3 2 d(3 x d/3 2 



u — xb 

XB 



5/2 



v 1+3 ^ 2 e( Wr - w) K - w) e/2 w e/2 



(18) 



The angles (3\ and (3% are the polar and azimuthal angles of one of the unobserved partons defined in the mentioned 
frame. The orientation of the axes in this frame was chosen in order to simplify the functions to be integrated, v and 
w are the energy and the angle of the fragmenting parton, respectively. In this case, as we anticipated, there is no 
correlation between them, but the 9 function splits the integration region R in the u, v and w volume: R — B0UB1UB2 
where BO is given by 



BO = {u e [xb,%u], v £ [vh, a], w G [0, 1]} . 



(19) 



and Bl and B2 are the LO regions given in eq.©- 

As it is common practice, in order to perform the angular integration, matrix elements have to be decomposed, via 
partial fractioning, in such a way that all the angular integrals end in the standard form |T^ 



I(k,l) 



dPi / dp 2 



sm 



X+e 



/?i sin e 2 



(a + bcos/3 1 ) k (A + B cos fa +C sm fa cos (3 2 ) L 



(20) 
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This kind of integrals can be classified into four categories according to whether their parameters satisfy either a 2 = b 2 
or A 2 = B 2 + C 2 , both relations simultaneously, or neither of them. In the present case, after the partial fractioning 
we obtained 31 independent integrals. 23 of these integrals were calculated to all orders in e extending the results of 
refs. 0,0,^3- The remaining 8, which we were not able to calculate to all orders, need to be carefully handled 
before expanding them in a power series in e. 

The difficulty with the above mentioned integrals is that e is not only regulating the (3 integration, but also the 
singularities in the remaining variables: u, v and w. Although the integrals may be regular functions of these variables, 
their coefficients may be not. An illustrative example of this situation is given by the integral 1(1, 1)\a 2 =b 2 +c 2 i as in 
eq. (I20|l with k = 1 , I = 1 , and satisfying A 2 = B 2 + C 2 , for which the order by order computation in e gives (eq. C30 



j(i,i) U2=B2+c2 = -^^{H + io g 



(aA-bB) 2 
[a 2 - b 2 )A 2 



0(e) 



(21) 



Let us first consider the case when a + b ~ (1 — w). As long as the coefficient of this integral is regular at w = 1, 
the above expression is integrable as function of w. However, if the coefficient has an extra factor of (1 — w)^ 1 the 
resulting expression is ill defined due to the logarithm in eq. l|21|l which behaves as log(l — w). In order to skip this 
problem, one can recast ea. l|20|l using the general methods described in the appendix A of ref.^l| obtaining: 



J(M) 



--B 2 +C 2 




A 



A(a-b) + b(A- 
B) - 2Ab 



B)x 



b{A-B)x 



(22) 



The integral in last equation can be splitted into two pieces, one containing the hypergeometric function, which can 
be integrated order by order in a power expansion in e after factoring out (a + by/ 2 , and the other which can be 
integrated to all orders in e: 



I(1A)\aw 



+c 2 



,/ i - hu i 7 + 21og 



;A-bB 



(a - b)A 



+ (a + b) €/2 



-- + bg[o-6] 



0(e) 



(23) 



Notice that factoring out (a + 6) c / 2 before the expansion in powers of e avoids the appearance of powers of log(a + b) 
in the series, which would be singular in the w — > 1 limit. In this way we obtain well defined integrals in w, as long as 
e > 0, even if the coefficient has a pole in w — 1, which is rather frequent. In some cases, the angular integrals have 
singularities in u, v or w by themselves, but they can be managed in the same way as in the example above. 

The procedure just illustrated was performed for all the 8 integrals and for the different combinations of singularities 
in u, v and w, expressions for them are available upon request. It is also important to stress that, as the singular 
distributions that show up in the matrix elements after the angular integration give rise to additional poles in e, it 
is necessary to calculate contributions up to order e 3 in the angular integrals. Fortunately, these poles are always 
accompanied by one or more 5 functions and those higher order terms only need to be calculated in the corresponding 
limits, what simplifies considerably the integrals. 








FIG. 5: One loop contributions to the a s gq cross section. Diagrams obtained from the first four by reversing the quark line 
must be also taken into account. 



Virtual contributions for the gq subprocess are obtained from the interference of the one loop graphs in figElwith 
the box graphs in figlSl Integration over the loop momentum was done using the standard Passarino-Veltman [l9| 
reduction algorithm and computing the resulting 2, 3 and 4-point scalar integrals. The integration over the phase 
space of the unobserved antiquark can be trivially performed using the energy-momentum conservation 8 function. 
After this integration, the remaining phase space can be written as 



dPS {2) = 



1 



87rr(l + e/2) 



91 

4/T 



72 



u (1 - x B ) f 1- Xb 



Xb 



XB 



s/2 



v e w e/2 (l - w) e/2 S(w r - w) dv dw , 



(24) 
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where v and w are the energy and angular variables of the hadronizing quark respectively. As can be seen from the 5 
function in ea. Q24[) . these two variables are correlated, which is a distinctive feature of the two particle phase space. 
It implies an additional constraint over the integration region: as w — w r < 1 then v > a should be satisfied and, as 
it happens for the order a s corrections, the integration region V has to be splitted into V = Bl U B2 where Bl and 
B2 have been already defined in eq.l©- 



III. SINGULARITY STRUCTURE 



Once the angular integrations are performed, the hadronic tensor shows a rich variety of singularities in the (u, v, w) 
space, regulated by the parameter e. As it is standard in this kind of calculations, the above mentioned singularities 
should be prescribed in order to get a series expansion in powers of e suitable for making explicit their cancellation. 
These cancellations are performed by coupling constant renormalization for the UV singularities, by cancellations 
between virtual and real contributions for the soft ones, and by renormalization of parton densities, fragmentation 
and fracture functions in the collinear case. 

A standard example for the above mentioned prescriptions is the appearance of factors like (1 — u)~ 1+£ in the totally 
inclusive cross section where, after the phase space integration, u is the only remaining variable. In this case one can 
use the standard substitution: 

(l-u)- 1+ ^ = -S(l-u)+(-^—) +0(e), (25) 



e V 1 — u 

where (1/(1 — u)) +u <p ii is the usual 'plus' distribution 



h«[o,i] 



V_!J) m .r du m^m. (26) 



o \1 u /+u[0,l] Jo I u 

However, in the one-particle inclusive case, the structure of the singularities is much more complex, mixing the three 
variables and consequently this simple prescription is no longer adequate. In fig.© we show the curves along which 




FIG. 6: Position of the singularities in the v-w plane for xb < u < x u . The bold lines represent the curves where the hadronic 
tensor becomes singular. 

the singularities in the regions BO and Bl appear in the v-w plane after the angular integration is performed. We will 
focus on this two regions because they contain all the singularities in the forward direction which need to be factorized 
in the redefinition of fracture functions. The case of region B2 is quite similar to that of Bl without the complications 
of the poles in w = w r = 1 but with additional singularities along the plane u = 1, and it will be discussed at the end 
of this section. 
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A simple inspection of ng.JBJ allows one to distinguish different possibilities for the singularity structure of the 
terms in the hadronic tensor. In principle the integration leads to terms that can have none, one, or two poles along 
the thick curves in the figure, respectively. 

The case of a single pole can easily be handled with minor modifications to the prescription formula in ea. (|25[l . For 
terms with more than a single pole, the singular curves can either intersect themselves (for example poles in w = 1 
and w — w r ) or not (for instance w — and w = 1). 

Terms with poles along two non-intersecting curves can be shown to be always transformed by partial fractioning 
into two terms with single poles, which reduce to the previous case. On the other hand, the case of two intersecting 
singular curves can not be reduced to a simpler one and needs to be treated in a more subtle way. 

Overlapping singularities as those mentioned in the previous paragraph can be further classified according to 
whether: a) both curves lie in the integration region, like w — and w = w r in region Bl; b) one of them comes from 
the outside of the integration region but intersects it at some point, as it is the case of w = 1 and w — w r in BO; and 
c) both curves converge into a single point of the integration region but coming from the outside, like w = — (1 — v)/v 
and w = —xb(1 ~ v)/v. 

The first and third occurrences can be cast into the second, by partial fractioning, leaving us with only one case. The 
technique we employed to treat it is better illustrated by means of an example. Let us consider the two dimensional 
integral 



1(e) = [ dx f 1 dy f(x, y) (1 - y)- 1+e (1 - x , (27) 

Jo Jo 

where f(x,y) is a regular function in all the integration region. The integrand has poles along the curves y = 1 and 
y = 1/x which intersect at x = 1, y = 1. These singularities are regulated if e > (notice that the integral remains 
finite even if the term e in the exponent of (1 — xy) is absent). If one wrongly uses the recipe in ea. l|25(l to prescribe the 
singularity in y = 1 and then again to deal with the pole in x = 1 coming from the 5 term, one ends with ill-defined 
terms (more precisely terms with 'plus' distributions which are not integrablc) and the leading singularity, in this case 
a double pole e~ 2 , is accounted twice. The correct way to deal with this integral is to re-write it as 

1(e) = f dx f 1 dyf(x,l)(l-y)- 1 +<(l-xy)- 1+f - 
Jo Jo 

r dx f dy (f(x,y)-f(x,l)) (l-y)- l+ '(l-xy)- 1+ *. (28) 
o Jo 

The second term is integrable in the limit e — > whereas in the first one the integration over y can be performed and 
gives 

He) = I 1 dx (1 - x )-»* f(x, 1) ^ 2 M + ^] + f 1 dx r dy ^ y \- f{xA ] + O(e) . (29) 



o 



(1 - y)(l -xy) 



Now, the integral in the first term can be prescribed using Q25[l. Doing that substitution, we end with the following 
identity: 

}S(l-y) 




1 



(i - y)0- ~ x v) 



[0,1] J 

0(e) . (30) 



y[o,i] 



The 'plus' distribution 1/((1 — y)(l — x y)) y [o,i\ stands for the second term in the r.h.s. of ea. l|29(l . The factor 1/2 in 
the double pole is a consequence of the fact that the singular curve y — 1/x only intersects the integration region in 
a single point. 

Prescriptions for all the singular (but regular at u = 1) terms appearing in the matrix elements can be found, 
besides some subtleties related to the integration intervals in Bl and B2, with the technique shown in the example. 
Expressions for the prescriptions relevant in the w = 1 region can be found in Appendix II. 

The only remaining item is the prescription of singularities in u = 1 in B2. These poles always appear as factors 
1/(1 — u) and only give rise to singular integrals (when e — > 0) in terms proportional to 5(w) or 5(w r — w) that come 
from the prescription of the singularities in the v-w plane. This is so because of the upper limit w r in the w integration 
which goes to zero when u — > 1. For the 5 terms, the prescription of the singularities in u = 1 can be done exactly as 
in ea. l(25|) . 
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IV. FACTORIZATION OF SINGULARITIES 



As we mentioned in the previous section, once the angular integration and the prescription of the singularities in the 
u, v and w variables are accomplished, the partonic cross sections exhibit a complex structure of poles in e. Explicit 
expressions for this structure in region BO can be found in Appendix III. Adding virtual and real contributions all 
IR divergences cancel out, leaving us with the UV and collinear poles. UV poles are canceled by means of coupling 
constant renormalization: 



a. = a s {Ml) ( a s (M R ) fa ( M% ^ 



/r- 73 . (31) 



where Mr is the renormalization scale and /3q is the lowest order coefficient function in the QCD (3 function: 

0o = y Ca ~ t UfTp (32) 

with Ca—N for SU(N) and Tp = 1/2 as usual; nj stands for the number of active quark flavours. 

Collinear singularities have to be factorized in the redefinition of parton densities, fragmentation and fracture 
functions. The redefinition of parton densities is exactly the same as in totally inclusive DIS whereas fragmentation 
functions are renormalized as they are in one-particle inclusive electron-positron annihilation. Expressions for renor- 
malized parton densities and fragmentation functions, up to order a 2 and in the MS factorization scheme, can be 
found in refs.0 and |2jj respectively. 

Notice that the renormalization of parton densities and fragmentation functions implies convolutions between the 
evolution kernels and the SIDIS cross sections. At variance with the totally inclusive case, the convolutions between the 
0(a s ) cross section and the LO kernels include plus distributions in more than one variable which need to be handled 
with care. In order to make explicit the cancellations between the 0(a 2 ) cross sections and these counterterms, the 
results of the above mentioned convolutions need to be expressed in terms of the very same variables used for the cross 
sections. One way to accomplish this is to retain to all orders in e the 0(a s ) cross sections, that is without replacing 
the singular factors like (1 — in terms of distributions as described in the previous section, and rewrite the 

plus distributions in the LO kernels using 



lim (1 - x)- 1+e ' --5(1- x) + O(e') . (33) 



1 ~ X J+x[0,l\ e '^° 

In this way the appearance of plus distributions is avoided and the convolutions can be explicitly performed. The 
resulting expressions can be prescribed, keeping up to constant terms in e and e', in exactly the same way as the 
0(a 2 ) cross sections. Notice that at this point the poles in e' must cancel and the limit e' — > can be safely taken, 
reflecting the fact that the LO kernels were already regular. The above mentioned procedure allows to extend the 
results of ref.|13 to SIDIS. 

Once the renormalization of parton densities and fragmentation functions is accomplished, the remaining singulari- 
ties occur in the region BO and are proportional to 8(1 — w), that is the forward direction, so they have to be factorized 
into renormalized fracture functions. Otherwise, factorization would be broken. The bare fracture functions can be 
written in terms of renormalized quantities as: 



-jfaMriMj^pfi^Mf) (34) 



s u \u 

where the factorization scale has been chosen to be the same for the three distributions. The functions Aj»_j and 
A/.j«_j are fixed in order to cancel all the remaining singularities in the cross section. 

The homogeneous kernels Aj<_j are the same that appear in the inclusive case for parton densities and can be 
obtained from the corresponding transition functions in ref. |10| . whereas the non- homogeneous A/-j<_j are presented, 
for the case j = g, in this paper for the first time. Explicitly: 

/ o \ e/2 

= ~fr U ~P^ 9 M (35) 
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(u,v) = A g ^ g (u,v) = (^)\f 2 (-^) ||(pi,l,(«,i;)®Pi2 ff («) + ^L fl («,«)(»pW,(«) 



(36) 



M?\ £/2 



2 ^ 



47T^ 2 



47T^ 2 



-o (Pil„(u,t;)® A («) 



(37) 



where a s is the bare coupling constant, the convolutions are defined as 

f(u,v)®g(u) = [ ^f(u,v)g(V), 

f " dv _ ( v \ 
f(u,v)®g{v) = J —f{u,v)gy-j , 

" 1 ~ uv du 



(38) 



f(u,v) ®' g(u) 



du u ( _ u \ /m\ 
— - J [u, -v \g - 
u u \ u ) \u) 



and 



P£l j (u,v) = P£l j (u)Hv- 



_ p(0) 



1 - u 



(39) 



Finally the 0(a 2 ) kernels are given by 



P£La(v,v) = C A T F \ - 



(3-8u)u 4(l-u)u 



8 m 3 



!m 2 (1 + m) 2m(1 + 4m-3m 2 ) 



(1 — MI)) 4 (1 — «I)f 



(1 — u v) 



2(1-3m)m , 
^ — + log 



+ log 



1 — uv 

1 — u — u v 



v \ 2 P^ g (u) 



v 



1 — u J v (1 — uv) 

6P 9 ( 2 9 (m) 



+ log(l + u) 



„ \ 2 2mR ( 2„(-m) 
-u (l + 2u)-u 2 v + 



1 — MI) 



2 (1 - 3u) m - 



„ 2 2m (1 + 4m)2m (1 + 2m + 4m 2 ) 

-3m 2 m-+- 



+ 



4 m 2 (1 + u) 4m 3 



+ log(u) 



(1 — uv) 3 (1 — iiu) 4 

9 p(°) / ^ 
u (3 + 2m) + " W + m 2 m + 



4 + m + 2m 2 v + 
8 m 3 



1-U!) 



, , 4P<2 3 (m) (2-uv) 



(1-un) 

+ log(l — uv) 



v (1 — uv) 

u 2 (l + u) 4 m (1 + 2m + 4m 2 ) 4m(1 + 3m) 



(1 — uv) 4 (1 — uv) 3 



(1 — Ml))' 



+ 



C F T F |4u 

m 2 (2 + v) 
1 - m 

2m 4 m 



m 2 (l-4u) u 3 v 3 



2 + 5m 

1-u (1-w)" ' (l + vj 2 T+v 



+ log(u) 



+ 



1-wti 
1 1 - 2u 



+ log(l + m) 



(i-«r 



+ 4P<2,(«) 



4m 3 



1-u f (1-u)" 



2 (1 - 2u) „ 2 ,„ 

— i- '- + 2u 2 (2 + v 

1 + v 

\- log(l — u — uv) —2m — 



(1-u) 2 (1 + v) 2 1 + v 
-4 m 3 



+ log(l - u) 
2m 3 w 



1 -m 
1 



(1-M,) 2 (1+M) 2 
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1 + v 



1-M 



and 



P^lLi^v) = C A T F { -m(1-4m)- 



V 1-U (1-M) 2 y 

8(1 + m) 2(1 + 4m-3m 2 ) 2(1-3m) 



(40) 



M (1 + I)) 4 m(1+w) 3 



m(1 + m) 2 



+ 



1 + v 



+ log(l — uv) 
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(l-2uv + u 2 {l + v 2 )) ,P!j d 9 (u) , , 

— — + log(l - u) , + log 



1 + V 



1 + v 



1 — u — U V 



3 u 



log(u) 



„ 2 2wP„ ( 2„(w) 
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log 
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u (1 + w) 4 u (1 + w) 3 
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(i-uY 
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where at = (1 — tt)/ti. Although A gg ^ 9 is formally a NLO kernel, it occurs for the first time at order a 3 thus it does 
not show up in the present calculation. Notice that the NLO kernels depend on both it and v variables and that this 
dependence cannot be factorized. 

Once obtained the explicit expressions for the relations between renormalized and bare fracture functions, we can 
easily derive the evolution equations for the renormalized fracture functions, which can be written as 



dMl h/p (t[,(,M 2 ) _ as{ M 2 ) f 1 du 



d log M 2 

a s {M 2 ) 1 fiTZ du 

u 1< 



2n 



2tt £ 



dv 

v 



e u 



P£i J (u,v) + ^-^P^(u,v) 



M. 



as(M 2 ) p(1) 
2tt 



c 



where the NLO kernels P^-(u) arc 1/8 of those given in ref.|22j due to the different conventions implemented. At 
variance with the LO case where the kernels are proportional to S(v — (1 — u)/u), the NLO kernels have support in 
all the integration region in the non homogeneous term of ea. (|42|l . Due to this fact, at NLO, the non homogeneous 
terms in the evolution equations do not take the familiar form given in eq.(12) of ref.[lj. In terms of moments, eq.(@2J) 
can be written as: 



d log 

where the moments are defined as 



= Kh/P K n] Pi^j [m) + fr /P [m + n- 1] D r h/k [n] P kl ^ [m - 1, n] 



F[m,n] = fjf jTCF^O, F[m] = jf* | F(£) , 



(43) 



(44) 
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f = 0. 1 f=0.1 ( = 0.0' f=0.1 (=0.001 ( = 0.1 




log(QVGeV') log(Q 2 /GeV) log(QVGeV') 

FIG. 7: Non homogeneous contributions to the derivative of M q for different values of £ and f. Inset plots show the integral 
over Q 2 of this contributions taking Mq' H ' (Qo) = with Qo = 1 GeV as a reference. 

and 

Pu^j [m, n] = J 1 ^ J 1 6 ^ r C" iWj (e, (45) 

Figure [7| compares (for different values of £ and £) the relative size of the LO and NLO contributions to the 
non homogeneous term in the evolution equation Ij42(l computed with standard sets of parton distributions |2ia] and 
fragmentation functions [24j for the case i = q and h = n + . The inset plots show the integral over Q 2 of this 
contributions. Notice that only those terms proportional to f g /p were taken into account in the 0(a 2 ) pieces. In 
reference || it was found that the LO non homogeneous contribution falls rapidly as £ grows. This behavior is related 
to the shrinkage of the integration region and with the fall of fragmentation functions, D^/i^z), in the limit z — ► 1. 
This is also the case of the NLO contributions. At moderate and large values of £ (£ > 0.1) the 0{a 2 s ) contributions 
are typically one order of magnitude smaller than the 0{a s ) ones so NLO and LO results differ only by a few percents. 
This can be traced back to the extra power of a s and the small size of the integration region since the interval of 
the v integral in ea. (|42|l shrinks to the point (1 — u)/u when £ — ► 1 — £. However, when £ diminishes the integration 
region expands and NLO contributions grow considerably faster than the LO ones which are kinematically restricted 
to the curve v = (1 — u)/u. The remarkable growth of the 0(a 2 ) terms makes these contributions even larger than 
the constrained 0(a s ) pieces at lower values of £ and thus a priori non negligible in the evolution equations. 

Of course, in order to assess the actual relevance of the NLO non homogeneous effects in the full evolution of 
fracture functions, one needs a realistic (based on actual data) estimate for the size and shape for these functions 
at a given scale, and compute the evolution taking into account all the appropriate kernels, but our present results 
suggest that non homogeneous NLO effects could be relevant. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have computed the 0(a 2 ) gluon initiated QCD corrections to one particle inclusive deep inelastic 
processes. At variance with the inclusive case, in one particle inclusive processes the kinematical characterization 



13 



of the final state particle requires to preserve the full dependence of the amplitude in the relevant variables. This 
impedes the cancellation of some singularities to be later factorized into fracture functions and leads to a non trivial 
singularity structure. In order to deal with this we have highlighted the importance of collecting to all orders the 
potentially singular factors in the 3-particle final state angular integrals and implemented a general approach for the 
prescription of overlapping singularities. 

By the explicit replacement of the bare parton densities, fragmentation and fracture functions with the corresponding 
renormalized quantities in both the 0{a s ) and 0(a 2 s ) cross sections, we have explicitly verified the factorization of 
collincar singularities obtaining for the first time the relevant kernels at this order. In doing so, we give a recipe 
for dealing with convolutions of distributions in more than one variable which occur in the computation of the a 2 . 
contribution coming from the convolution of 0(a s ) cross sections and renormalized functions. We also derived the 
evolution equations for fracture functions valid at NLO. 

Regarding the phenomenological consequences of these corrections, we have found that the 0(a s ) contributions to 
the evolution equations are mild in most of the kinematical range, however they are as important or even larger than 
the a s ones for small values of xb, where these last contributions are suppressed by the available phase space. This 
behaviour, at variance with the LO case, allows the non homogeneous effects to be sizeable even at larger hadron 
momentum fractions, thus being relevant for the scale dependence for diffractive and leading baryon processes. 
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Appendix I 

In this appendix we present the results obtained for the coefficients in the 0(a s ) cross sections. They are 
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Appendix II 

The redefinition of fracture functions can only factorize singularities in the forward region. In the hadronic tensor 
these singularities show up after the angular integration as (1 — w)~ 1+£ factors which have to be prescribed as 
explained in section III. As we mentioned there, special care has to be taken with overlapping singularities. In this 
case, inspection of figure [5] shows that the only problematic configurations are terms singular along w = 1 and w = w r 
in regions BO and Bl. Using the procedure described in section III we obtained suitable prescriptions in both regions: 
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where ei and e 2 are multiples of the regulator e. Notice that these expressions are valid to all orders in e. Terms 
singular only along w = 1 can be prescribed using the standard rule in eq. I|25[) . Terms singular in w — 1 and w = 
can be managed by partial fractioning: 
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and be prescribed also as in eq. (J^SJ. 



Appendix III 

The singular pieces of the order oq partonic cross sections a gg and b gq in region BO can be written as 
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where the Pj^j_j{u, v) are defined in ea. (|39[) and the functions C\ ^(u) are the coefficient functions of totally inclusive 
DIS at 0(a s ), they can be found in refs. [aliof. Convolutions between kernels are as in eq. I|38|l with the replacement 
v —*■ (1 — xb)v/xb whereas the convolutions between kernels and coefficient functions are given by 
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